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ABSTRACT 

Mass is the most important stellar parameter, but it is not directly observable 
for a single star. Spherical model stellar atmospheres are explicitly characterized 
by their luminosity (L^,), mass (M*) and radius (-R*), and observations can now 
determine directly L^, and i?*. We computed spherical model atmospheres for 
red giants and for red supergiants holding and i?* constant at characteristic 
values for each type of star but varying M^, and we searched the predicted flux 
spectra and surface-brightness distributions for features that changed with mass. 
For both stellar classes we found similar signatures of the star's mass in both 
the surface-brightness distribution and the flux spectrum. The spectral features 
have been use previously to determine logiQ^g), and now that the luminosity and 
radius of a non-binary red giant or red supergiant can be observed, spherical 
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model stellar atmospheres can be used to determine the star's mass from cur- 
rently achievable spectroscopy. The surface-brightness variations with mass are 
slightly smaller than can be resolved by current stellar imaging, but they offer 
the advantage of being less sensitive to the detailed chemical composition of the 
atmosphere. 

Subject headings: stars: atmospheres - stars:fundamental parameters - starsdate- 
type - stars: supergiants 



1. Introduction 

Plane-parallel model stellar atmospheres are characterized by their effective tempera- 
ture, defined as 

and surface gravity. 

As is apparent, these two parameters are degenerate expressions of the three more funda- 
mental stellar parameters L^,, and i?*. Along the main sequence the variation of i?* is 
small, making Teg and log^glS') nearly unique identifiers of the star. However, in the cool, 
luminous quadrant of the Hertzsprung-Russell diagram, stars with a large range of i?^ and 
converge to similar values of Tgfj and logiQ^g), which makes it difficult to establish the 
fundamental parameters of non-binary red giants and red supergiants. 

Fortunately, optical/infrared interferometry is helping to break this degeneracy by mea- 
suring directly the angular diameters of stars. By combining these angular diameters with 
stellar parallaxes that are being measured with increasing precision to ever greater distances, 
it is now possible to determine i?* observationally. In addition, combining the stellar flux 
measured at Earth with the same stellar angular diameter and parallax yields the L^,, as- 
suming the star radiates isotropically. With measured values of L^, and i?^, the observational 
attention turns to M*, the key parameter that determines the other two. 



It is now routine (IHauschildt et al.lll999l : iGustafsson et al.ll2008l : iLester fc Neilsonll2008l ) 



to compute spherically symmetric models of stellar atmospheres that are characterized ex- 
plicitly by L^, and i?^, in addition to the chemical composition. Assuming that red 
giants and supergiants are spherically symmetric, the ability to measure and Ri, raises 
the question whether spherical stellar atmospheres predict any indicators of stellar mass that 
could be determined observationally. 
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Red Giants 



To search for possible indicators of stellar mass, the SAtlas program (ILester fc Neilson 



20081 ) was used to compute a spherical model atmosphere with properties representative of 



a typical red giant: = 500 Lq, = 0.8 Mq and R^, = 50 Rq. The model was calculated 
with the opacity distribution function (ODF) version of SAtlas using the ODFs recomputed 
by Castelli (wwwuser.oat.ts.astro.it/castelli/) in 2012 February to incorporate Kurucz's cor- 
rections to his H2O line list. As Kurucz explains (kurucz.h arvard.edu/molecules/H20/h2ofastfix.readme), 
he uses the gf values from [Partridge fc Schwenkd (119971 ) after renormalizing them to agree 
with the gf values from the HITRAN (High-resolution TRANsmission) molecular absorp- 
tion database (www.cfa.harvard.edu/hitran/), but ignoring nuclear spin factors. Kurucz's 
original H2O line list did not take into account that the Hg^O, Hg^O, H2^0 and HD^^O iso- 
topologs each has its own partition function. To correct for this, while also using only the 
H2^0 partition function in Atlas, Kurucz assumes that, at stellar temperatures, the par- 
tition function of each isotopolog differs from the partition function of Hg^O by a constant 
multiplicative factor, which he det ermines by comparing the gf values of a few strong lines 
in the HITRAN database with the [Partridge fc Schwenkd (119971 ) values. The ODF assumes 
a depth- independent microtu rbulent velocity of 2 km/s and the s olar abundances tabulated 
by iGrevesse fc Sauvall (119981 ) . The mean differ ence between the iGrevesse fc Sauvall (119981 ) 
abundances and the more recent tabulation of lAsplund et al.l (120091 ) for the elements with 
number densities > 10~^° x Ah is only 0.02 in log^o of the abundances. The ODF also uses 
the large wavelength spacing with an average spectral resolving power of R = A/AA ~ 50, 
corresponding to AA ~ 10 nm at 500 nm. A test using an ODF with the smaller wavelength 
space, which gives ^ 5x higher spectral resolving power, produced a model structure (tem- 
perature as a function of gas pressure) that was almost indistinguishable from the structure 
of the model with the lower resolution ODF. Therefore, the lower resolution ODF was used 
to calculate all the models. The model was iterated until the luminosity was constant to 
< 1% for lo gio('^Ross) > 10 ^ and the derivative of the luminosity was constant to < 1% for 
log]^o(''"Ross) > 10~^, where tross is the Rosseland optical depth. 

To test the dependence on mass, two additional red giant models were computed holding 
the values of and R^ constant, but with = 1.6 Mq and = 2.4 Mq. These 
parameters correspond to Tes = 3865 K for all three stars, but log]^o(fi') = 0-94 for the 
0.8 Mq star, hg^oig) = 1.24 for the 1.6 Mq star and logio(5') = 1-42 for the 2.4 Mq star. 
The effect of mass on the stellar atmosphere can be represented by a relative extension 
parameter, e, defined as the radial distance from the deepest visible layer of the atmosphere 
out to an arbitrary upper radius, normalized by the stellar radius, 

_ i?(upper) — i?(deep) 



R. 



(3) 
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A convenient value for i?(deep) is where tross = 1, and the radiation's probabihty of escape 
is exp(— 1) = 0.36. The definition of the upper radius is arbitrary. Defining i?(upper) where 
tross = ~3, corresponding to about three pressure scale heights, leads to e(— 3). Using this, 
the 0.8 Mq star has e(-3) ^ 0.025, the 1.6 Mq star has e(-3) ^ 0.012 and the 2.4 M© 
star has e(— 3) ~ 0.008. Defining the upper R is where tross = —6, leads to the definition 
of e(-6). Using this, the 0.8 M© star has e(-6) f« 0.06, the 1.6 M© star has e(-6) ^ 0.028 
and the 2.4 Mq star has e(— 6) ~ 0.019. Clearly the extension parameter is greater with this 
choice of i?(upper), but by either definition these test models have modest extension. 



2.1. Surface Brightness Distribution 

The three red giant model atmospheres were used to calculate the surface intensity at 
100 values of the disk's fractional radius, r/i?^, ranging from to 1 in steps of 0.01, for each 
wavelength bin of the ODF that had a non-negligible amount of radiation. A comparison was 
made of the surface brightness distributions of the same model parameters computed with 
both the lower and higher spectral resolution ODFs described in the previous section. At 
wavelengths > 1000 nm, the surface brightness distributions were almost identical for the two 
spectral resolutions because there are so few spectral lines. At wavelengths < 1000 nm, the 
surface brightness distributions do differ very slightly because there are more and stronger 
spectral lines. The ODF with the finer wavelength spacing resolves some spectral detail that 
is averaged out in the lower resolution ODF. As a result, at some wavelengths the surface 
brightness distribution is shifted slightly toward or away from depending on whether the 
higher resolution ODF happens to include a strong line or not. The lower spectral resolution 
surface brightness distribution are used here because they seem to be a better match to the 
observations, but the synthetic spectra described in § l2.2l can compute the surface brightness 
distribution for any spectral resolving power, which can match any particular observation 
exactly. Figured] shows the limb darkening curves for the three red giant models at t he ODF 



wavelengths close to the effective wavelengths of the V and R magnitudes listed bv iBessell 



( 20051) and to the effective wavelengths of the /, J, H and K magnitudes given by Tokunaga 



fcooof i. [t is apparent that increasing the stellar mass makes the relative intensity greater 



at every location on the stellar disk for every wavelength. This trend results from the gas 
density and the thermal radiation of the atmosphere increasing with the stellar mass. It 
is also apparent that the curves flatten as the wavelength increases toward the peak of the 
Planck distribution for the star's temperature. 

The differences between the curves in Figure [1] uses a knowledge of the stellar limb's 
absolute location, i?^, which is present in the models. However, as the stellar mass is 
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Fig. 1. — Surface brightness distributions for three red giant spherical model atmospheres 
at several wavelengths. The atmospheres all have = 500 Lq and i?^ = 50 -R0. The solid 
lines are for = 0.8 Mq, the dashed lines are for = 1.6 Mq and the dotted lines are 
for = 2.4 Mq. The wavelengths are approximately the effective wavelengths of the V , 
R, I, J, K and H filter systems. 
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reduced, the upper atmosphere becomes more tenuous, and the radiation from these outer 
radii decreases to the point where some of the curves in Figure [1] appear to go to zero before 
reaching i?^. In fact, the relative intensity is not zero, but it is less than 0.01, making it too 
small to be resolved in the plot. This reduced relative intensity at the stellar limb is a direct 
consequence of the spherical geometry of the models, and it also affects the observations. 

To recast the differences in Figure [1] into a form that is suitable for observational ap- 
plication, we define the location of the stellar limb to be where the surface brightness has 
fallen to a specified fraction of the central brightness, 

I(r) 

-Riimb = r where — — - = /. (4) 

I[r = 0) 

Using this definition of Rumh, we renormalize the surface brightness distribution. Figure [2] 
shows the result of setting / = 0.01, 0.02, 0.05 and 0.1 for one wavelength of one spherical 
model. The wavelength A = 705 nm is approximately in the middle of the range shown in 
Figure [H From an examination of Figure |2] we elect to define -Rumb using / = 0.01. This 
choice gives a value of -Rumb that is close to while also being within reach of observations, 
but clearly another value of / could be chosen. Figure [U shows that the surface brightness 
falls off more steeply at the limb as the model's mass increases, which reduces the difference 
between -Rumb and R^,. With an angular resolution ^ 1% of the stellar disk and a photo- 
metric precision of ~ 1%, observations of the surface-brightness distribution from optical 
interferometry or a planetary transit could determine this location of the stellar limb. 

Using this observationally motivated definition of the stellar limb, we renormalized the 
intensity curves at each ODF wavelength and searched for the fractional radii where stellar 
mass produced the largest intensity difference. As an example. Figure [3] shows the differences 
between the r/i?iimb values for the red giant models with 0.8 Mq and 1.6 Mq and between 
1.6 Mq and 2.4 Mq models for the various limb renormalizations considered in Figure [21 The 
wavelength A = 705 nm is typical of a wide range of wavelengths. The difference between 
the 1.6 Mq and 2.4 Mq models is nearly the same for all the limb renormalizations, with the 
difference being greatest near /(r)/J(r = 0) ~ 0.5. The difference between the 0.8 Mq and 
1.6 Mq models is also greatest near I{r)/I{r = 0) ~ 0.5 for the renormalization at / = 0.01, 
but as the renormalization point is located farther in from i?^, the location of the greatest 
difference between the two models also moves inward. We elect to use the difference at 
/(r)//(r = 0) = 0.5 because it applies to the limb normalization with / = 0.01 and because 
the half-intensity point is easy to determine observationally. 

We determined the renormalized fractional stellar radius, r/i^umb, at the half- intensity 
point for each ODF wavelength from 500 nm to 2000 nm, which is shown in Figure HI The 
models show similar trends with wavelength. For < 900 nm the half-intensity point moves 
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Fig. 2. — Renormalization of the limb-darkening profile at A = 705 nm for the spherical 
model atmosphere with = 500 Lq, = 50 Rq and = 0.8 Mq. The different curves 
represent different choices for the value of "/" in Equation H] to define -Rumb- The solid line is 
for / = 0.01, the dotted line is for / = 0.02, the dashed line is for / = 0.05 and the dot-dash 
line is for / = 0.1. 




I(r)/I(r = 0) 



Fig. 3. — Difference between r/i^nmb at A = 705 nm plotted as a function of tlie renormalized 
surface-briglitness distribution. Tlie solid line is for the limb defined using Equation H] with 
/ = 0.01, the dotted line is for / = 0.02, the dashed line is for / = 0.05 and the dot-dash line 
is for / = 0.1. The top plot shows the difference between the 0.8 Mq and 1.6 Mq models, 
and the bottom plot shows the difference between the 1.6 Mq and 2.4 Mq models. 
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Fig. 4. — Top plot displays for the three red giant spherical models the wavelength depen- 
dence of the renormalized fractional radius (r/i^umb) where the surface brightness is half the 
central brightness. The spherical models all have = 500 Lq and i?* = 50 Rq. The solid 
line is for = 0.8 Mq, the dashed hne is for = 1.6 Mq and the dotted line is for 
= 2.4 Mq. Bottom plot is the differences between the models in the top panel, with the 
solid line being the difference between the 0.8 Mq and 1.6 Mq models and the dashed line 
is the difference between the 1.6 Mq and 2.4 Mq models. 
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toward the limb as the wavelength increases, but for A > 900 nm the fractional radius of 
the half-intensity point is nearly constant. Within these general trends the three masses are 
offset from each other, with the location of the half-intensity point shifting toward the limb 
as the mass increases. It is also clear from the bottom panel of Figure H] that the location 
of the half-intensity point varies significantly with wavelength for the lowest mass model 
considered here. For example, the separation is relatively larger at 705 nm, the wavelength 
used in Figure [21 than at adjacent wavelengths. The separation also increases near 1600 nm, 
where the continuous opacity due to H~ is minimum. 

Examining the stellar spectrum in the regions of larger separation, we find that they 
have relatively less line opacity. Because of this, the intensity at these wavelengths comes 
from the thermal emission deeper in the atmosphere. As the mass of the star increases, the 
greater gas density produces a larger amount of thermal emission. As a result, the surface- 
brightness distribution of the more massive star is flatter, and its half-intensity point shifts 
toward the limb. The conclusion is that it is important to measure the surface-brightness 
distribution with sufficient spectral resolution, i? ^ 50 — 100, to isolate the most sensitive 
spectral windows. Of the existing interferometers capable of working in the optical region 
where the separation between the curves is greatest, the CHARA array's VEGA instrument 



(iMourard et al.ll2009l ). which observes the band from 0.45 to 0.85 /im, has three choices 
of spectral resolution, R = 1700, 6000 and 30000, all of which are substantially greater 
than the required spec tral resolving power. The Navy Precision Optical Interferometer 



(lArmstrong et al.lll998l ) disperses the spectral band from 450 to 850 nm onto 32 channels. 



providing R ^ 50, just at the lower edge of the required spectral resolution. In the infrared 



the AMBER instrument on the VLTI flPetrov fc AMBER Consortium! 120031 ) could easily 



resolve the H~ feature. The bigger challenge for all the interferometers is achieving the 
angular resolution needed to locate the stellar disk's half-intensity point with an uncertainty 
of just a few percent. Although this appears to be beyond the current capability, the rapid 
progress being made in this active field makes the needed improvement seem likely. 



2.2. Spectrum 

Because the model's mass affects its structure, there can be associated variations of 
the model's spectral flux. To test this, we computed synthetic flux spectra covering the 
wavelength range from 400 to 1000 nm using the red giant spherical model atmospheres 
with My, = 0.8 Mq and 2.4 Mq. These spectra, which were computed using spherical 
radiative transfer with the spherical model structure, have a spectral resolving power of 
R = 10^ and include the atomic, ionic and molecular lines from the lists on Kurucz's web 
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site (1 http://k urucz . harvard . edu /linelists . html ) . Because we are comparing models to each 
other, not attempting to match observations, we have not included rotation or various kinds 
of broadening beyond the depth-independent microturbulent velocity of 2 km/s that is used 
in the ODFs to compute the model atmospheres. 

From a comparison of the synthetic spectra, we found that ~ 95% of the lines either 
remained nearly constant or became weaker as the stellar mass increased. For the largest 
weakening, which occurred in the 400-500 nm wavelength band, the line core decreased 
by about 0.15 in residual intensity for the mass range of our models. The amount of line 
weakening diminished steadily with increasing wavelength, with our models in the 900- 
1000 nm interval having line core residual depths that were smaller by no more than about 
0.07. 

At all wavelengths, the lines showing the greatest decrease with increasing mass were 
overwhelmingly ions of heavy elements, such as Nd II 536.1165 nm and La II 580.8313 nm. 
However, there is one interesting exception to this dominance by heavy ions, which is [O I] 
2p^ — )■ 2p^ with lines at 630.0304 nm and 636.3776 nm. Figure E] shows the beh avior 
of the 630.0304 nm line, which is consistent with the work of iBonnell &: Belli ( 1l993bl ) who 
found that [O I] is a useful indicator of surface gravity. Of course, log^gls') combines and 
R^, so it is a composite indicator, unlike our approach. For the 2.4 Mq star the equivalent 
width of [O I] 630.0304 nm is 23% less than for the 0.8 Mf^^ star , and the decrease of the [O I] 
636.3776 nm line is 31%. As shown by iBonnell fc Belli ( Il993bl ). comparing measurements of 
these [O I] lines with lines of the OH molecule near 1.6 fim yields the oxygen abundance and 
the surface gravity, or now the M^. 

Of the ~ 5% of the lines that became stronger with increasing mass, all of them were 
either MgH or SiH, with MgH strengthening by about twice as much as SiH. The largest 
increases in the residual line depth, about 0.08 for our models, were for lines of MgH in the 
X ^E"*" — > y4 ^H band head around 500 nm. An example of one line from this band is shown 
in Figure |6l where the line's equivalent width has become about 25% stronger going from 
the least to most massive model. 

The lines of MgH in the 500 nm spe c tral r egion have a long history o f serving as 
gravity indicat o rs, beg inning with lOhmanI ( 1l934l ). More recently, iBell et al.l (1l985l ) and 
Bonnell fc Belli (jl993al ) have explored the use of these lines to determine the gravities of cool 
giants. By using the MgH lines in combination with lines of Mg I, it is possible to determine 
first the Mg abundance and then the surface gravity. Again, our focus on mass instead of 
the composite surface gravity allows for a sharper distinction. 



As a test, we used the plane-parallel Atlas_ODF code ( iLester fc Neilsonll2008l ) and the 
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Fig. 5. — Rectified flux spectra of the line [O I] A = 630.0304 nm for the red giant spherical 
models having = 500 Lq and R^, = 50 Rq. The solid line is for — 0.8 Mq, the dashed 
line is for = 1.6 Mq and the dotted hne is for ^ 2.4 Mq. 
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Fig. 6. — Rectified fiux spectra of the MgH line at A = 500.3349 nm for the red giant 
spherical models having = 500 Lq and = 50 Rq. The solid line is for = 0.8 Mq, 
the dashed fine is for = 1.6 Mq and the dotted fine is for = 2.4 Mq. 
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same ODF file to compute models with values of Teg and logi^lg) corresponding to the red 
giant spherical models. We used these plane-parallel models to compute synthetic spectra 
with a spectral resolving power of 10^ using the same line lists employed above. The resulting 
equivalent widths of the [O I] and MgH lines, and their change with surface gravity, were 
nearly identical with values computed with the spherical model. This shows that these lines 
are sensitive to the mass, not the spherical geometry of the atmosphere. 

The other molecular lines in the spectrum, except SiH, show either almost no change 
or became weaker with increasing mass. To gain some perspective on this, we plot in the 
top panel of Figure [7] the number densities as a function of tross for several of the prominent 
diatomic molecules in the atmospheres of the red giant models with = 0.8 Mq and 
2.4 Mq. In the bottom panel of Figure [7] we show the variation of the temperature with 
tross ioY both models. The vertical dashed line in Figure [7] corresponds approximately to 
the atmospheric depth where the probability of radiation escaping to the observer becomes 
negligible. Note that only MgH and SiH have their peak number densities in the readily 
observed layers. This is a consequence of these two molecules having lower dissociation 
energies than an y of the other common molecules, and MgH, with a dissociation energy of 



just 1.28517 eV (jShayesteh et al.l 120071 ). is more than a factor of two less than any of the 
other molecules. 



2.3. Composition 

The spectral results discussed in § 12.21 clearly dep end sensitively o r i the c omposition of 
the star. It is for this reason that it was necessary for iBonnell fc Belli Jl993bh to use other 



spectral features to determine elemental composition before using either the [O I] or MgH 
lines to derive surface gravity. The surface-brightness distribution examined in § 12.11 also 
depends on the metallicity of the star, which is shown in the left panel of Figure [H As the 
metallicity is decreased from solar, the surface-brightness distribution becomes fainter at the 
stellar limb and brighter at r/R^^ < 0.9. Some of the relative intensities in the left panel of 
Figure [8] appear to go to zero before r reaches R^,. As explained earlier in the discussion of 
Figure [H with decreasing stellar mass the upper atmosphere becomes increasingly tenuous, 
causing I(r) /l(r=0) to drop below 0.01 as r approaches i?^. Even though the relative intensity 
is greater than zero, it is too small to be resolved in the plot. 

The change at the limb is obvious for the models, for which R^ is known, but the 
location of i?^ is not available observationally. Therefore, the same observational definition 
of the stellar limb introduced in § 12.11 is employed here, with the result shown in the right 
panel of Figure [3 There is now a clear progression of the surface-brightness distribution 



Fig. 7. — Top five panels show the variation of the number density of prominent diatomic 
molecules with atmospheric depth, represented by trqss, for two red giant spherical atmo- 
spheres, both having = 500 Lq and = 50 Rq. The solid lines are for Af^ = 0.8 Mq 
and the dotted lines are for M^, = 2.4 Mq. The vertical dashed line is at logiQ(TRoss) = —0.5, 
which is the approximate limit of our view into the stellar atmosphere. Bottom panel shows 
the variation of the temperature with optical depth for the two models, which provides 
additional context for the depth dependence of the molecular number densities. 
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Fig. 8. — The left panel has the surface-brightness distributions plotted as a function of the 
fractional disk radius, r/i?^, for three red giant spherical models, all having the same L^,, i?^ 
and M^. The solid line is the model with solar metallicity, the dashed line is the model with 
0.1 solar metallicity, and the dotted line is the model with 0.01 solar metallicity. The right 
panel plots the same surface-brightness distributions, now as a function of the renormalized 
fractional radius, r/i?iimb, defined by equation in § 12.11 The lines have the same meaning 
as in the left panel. 
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with metallicity, with the intensity curve steepening as the composition becomes more metal 
poor. Therefore, to determine the star's mass from its intensity distribution it is essential to 
know the star's composition. However, this requirement is much less stringent than what is 
needed to utilize a particular spectral feature, such as those highlighted in § 12.21 Instead of 
requiring a precise value of the composition for the element producing a particular spectral 
line, only the overall composition of the star is needed, which can be determined by analyzing 
a selection of spectral features, or perhaps adequately estimated from spectral classification. 



3. Red Supergiants 

To extend our exploration of mass indicators, we next considered red supergiants. Again, 
we use the SAtlas_ODF program to compute atmospheric models that are spherically 
symmetric, this time with the representative parameters = 50, 000 L©, i?^ = 500 Rq and 
solar composition. For the masses we chose = 8 Mq, = 16 Mq and = 24 Mq. 
Our final models achieved constancy in the flux and the flux derivative very similar to the 
red giant models. The fundamental parameters of these models yield T^s = 3865 K, the same 
as the red giant models, but log^Q^g) = —0.06 for the 8 Mq model, logiQ^g) = 0.24 for 16 Mq 
model and logio(5') = 0.42 for the 24 Mq model. Using the relative extension parameter 
defined in equation ©, the 8 Mq model has e(-3) ^ 0.027 and e(-6) ^ 0.065, the 16 Mq 
model has e(-3) ^ 0.013 and e(-6) ^ 0.031 and the 24 Mq model has e(-3) ^ 0.008 and 
e(— 6) ~ 0.020, which, again, are all modest extensions. 



3.1. Surface Brightness Distribution 

The surface-brightness distributions of the red supergiants were computed and ana- 
lyzed in exactly the same way as for the red giants, and the wavelength variation of the 
half-intensity point determined from the renormalized intensities is plotted in Figure HI 
Comparing Figure [9] with Figure H] it is obvious that the trends and features are very simi- 
lar, with just minor difference in detail. As before, the different masses change the surface 
brightness by a few percent, and for A < 900 nm a spectral resolving power of -R ~ 50 — 100 is 
essential to isolate these regions. For A > 900 nm the separation between the curves is more 
constant, and a lower spectral resolving power would be adequate to detect this difference. 
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Fig. 9. — Top plot shows for the three red supergiant spherical models the wavelength 
dependence of the renormalized fractional radius where the surface brightness is half the 
central brightness. The models all have = 50, 000 Lq, R^, = 500 Rq and solar composition. 
The solid line is the 8 model, the dashed line is the 16 Mq model and the dotted line is 
the model with = 24 Mq. Bottom plot is the difference between the models in the top 
plot, with the solid line being the difference between the 8 Mq and 16 Mq models and the 
dashed line is the difference between the 16 Mq and 24 Mq models. 
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3.2. Spectrum 

To search for signatures of mass, we used the spherical red supergiant model atmospheres 
with Mi, = 8 Mq and = 24 Mq to synthesize the flux spectrum from 400 to 1000 nm 
at i? = 10^ using the same line lists as for the red giants. Comparing the line depths, we 
again found that ~ 95% of the lines were either unchanged or became weaker. As was found 
for the red giants, the lines displaying the greatest weakening with increasing mass were 
located at the shortest wavelengths we computed. The trend found for the red giants, that 
the amount of weakening decreases with increasing wavelength, also holds true for the red 
supergiants, but the drop off is somewhat more gradual until the very longest wavelengths. 
The total variation over the spectral range 400 — 1000 nm is the same. 

As was found for the red giants, the lines that weaken the most with increasing mass 
were overwhelmingly ions of heavy elements, with the [O I] lines remaining noteworthy 
exceptions. The [O I] lines in the supergiants are about twice as strong as for the giants, 
but their percentage weakening with increasing mass is somewhat less. In addition, the 
equivalent widths of [O I] 636.3776 nm and 630.0304 nm for the red supergiants each weaken 
by comparable amounts, although the percentage changes are different, but for the red 
giants the equivalent width change of the 630.0304 line is almost 50% greater than for the 
636.3776 nm line. 

The < 5% of the lines that become stronger with increasing mass are all lines of either 
MgH or SiH, with the MgH lines changing by about twice as much as the SiH lines. Although 
the MgH lines of the supergiants are not as deep as they are in the giants, their depths 
change by about 20% more than between the most and least massive giant models, and 
this is also true for the change of their equivalent widths. The variation with depth of the 
number density of the diatomic molecules for the supergiant atmospheres is very similar to 
the variation with depth for the red giants shown in Figure [71 

4. Discussion 

Both the distribution of the surface brightness and the flux spectrum contain indicators 
of the mass of the red giants and supergiants. Of these, the spectral features of [O I] and MgH 
have been recognized and used previously, although they were associated with logio(fi'), which 
depends on the combination of two of the fundamental stellar parameters. With the ability 
to measure independently both and i?^, these lines can now be used to determine the 
mass of a non-binary star. Current spectroscopic capability is clearly sufficient to measure 
these lines with the precision required to do the analysis. 
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The major complication in using spectral features to determine the stellar mass is the line 
variability exhibited by some giants and by man y supergiants. In an ext ended observational 
study of Betelgeuse (a Ori, HR 2061, HD 39801) Icravl J200ol . l200ll . [2008h recorded variations 
in the strength of several neutral metal lines in a relatively small spectral window centered 
on A 625.0 nm. These variations were attributed to changes in the continuous opacity that 
alters the contrast between the line core and the adjacent continuum. Three-dimensional 
radiative hydrod ynamic modeling also p redicts time variation of the spectrum due to surface 
inhomogeneities (jChiavassa et al.l 120111 ). The spectral analysis presented here found that 
most lines, including the neutral lines studied by Gray, have a negligible variation with 
mass. The specific spectral lines that we identified probably do show the spectral variations 
caused by surface structure in addition to the dependence on mass. Therefore, to employ 
the mass-sensitive spectroscopic features found here, it clearly is necessary to normalize the 
line strengths to adjacent lines insensitive to mass in order to remove the variability caused 
by processes occurring within the atmosphere. This is in addition to eliminating the effect 
of the star ' s com p osition by normalizing to ot her lines of the s ame el ement, as was done by 



Bell et all ( llQSSf ). iBonnell fc Belli (Il993af ) and iBonnell fc Belli fll993bf l. 



The signature of mass in the surface-brightness distribution is just a few percent, which 
is smaller than the variation of the spectroscopic features. However, this small signature may 
not be as challenging as it appears from our analysis. In our tests the giant and supergiant 
masses were each varied by a factor of three, resulting in A logio(fi') — 0-48 in both cases, with 
the corresponding extension parameter e(— 3) < 0.03 and e(— 6) < —0.06 . Even careful, de- 
tailed observational analyses of red giants and supergiants frequently have uncertainties ~ 0.5 
or more in their determinations of \ogiQ{g). Arcturus provides an illustrative example be- 
cause it has data of the h i ghest quality, which have been analyzed repeatedly. The thorough 
study by iPeterson et al.l (119931 ) concluded that logiQ(5') = 1.5 ± 0.15, found by minimizing 
the differences between the observed high resolution flux spectrum and spectra synthesized 
from a grid of plane-parallel Atlas s t ellar a tmospheres. The equally thorough examination 
of Arcturus by I Griffin &: Lynas-Grayl (1l999l ) determined that logio(5') = 1.94 ±0.05 by min- 
imizing the differences between observed and computed spectral energy distribu tions, again 
using plane-parallel Atlas models. Recently, iRamirez fc AUende Prietol (|201l[ ) determined 
logio(5') ~ -'-■66 ± 0.05 for Arcturus using evolutionary isochrones. The agreement of the 
logiQ{g) values is just marginal considering the error bars. Stars that are less well observed 
than Arcturus often have significantly larger uncertainty in their \ogiQ{g) determinations. 
If the surface-brightness distribution could determine log]^Q(5') with an uncertainty of 0.5, 
it would be competitive with the currently employed methods. In addition, the physics re- 
sponsible for the surface-brightness signature of mass is particularly simple, just hydrostatic 
and thermal equilibrium. We did find that the intensity distribution depends on the general 
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metallicity of the stellar atmosphere, but estimates from spectral classification or a more 
quantitative determination from a selection of spectral lines would be adequate to determine 
this parameter. The resulting mass determination will provide a different perspective from 
the other methods that rely on significantly different and more complex physics. 

The measurement of stellar surface brightness distributions is a revolutionary obser- 
vational achi e vemen t, but one that is still developing. For example, the recent study by 
Lacour et al.l ( l2008l ) of Arcturus, a star whose disk is well resolved by optical/infrared in- 
terferometry, concluded that the disk was not uniformly bright, but was unable to observe 
exactl y how the surface brightness varied with position on the disk. Several studies of Betel- 



geuse flBuscher et al.lll990l : iTuthill et al.lll997 



Young et al.ll2000l : iHaubois et a 



detected large bright features on its surface, while other studies (lOhnaka et al 



2009[| have 
20091 12OI1I ) 



have found only slight deviation from a uniform disk, significantly smoother than the pre - 
dictions of the three-dimensional radiative hydrodynamic models (jChiavassa et al.ll2009l ). 
However, when observational agreement is reached about the true appearance of the sur- 
faces of these stars, a key point is that the features will be due to physical processes taking 
place within the atmosphere. These disk substructures will complicate the observations, but 
they do not invalidate the conclusion reached here that the basic trend of the stellar surface 
brightness distribution contains the signature of the fundamental global stellar parameter of 
mass. 



Neilson fc Lesterl (120111 . |2012| ) have recently explored a different, more statistical way of 
using the surface-brightness distribution of a cool, luminous star to determine its mass. Using 
the same S Atlas_ODF code, they computed thousands of spherical model atmospheres with 
L^, and i?^ thoroughly covering the parameter space of the red giants and supergiants. 
After computing the surface-brightness distribution for each model, a limb-darkening law 
with two coefficients was fit to the intensities. An analysis of the ensemble's limb-darkening 
coefficients showed that they depend on the atmospheric extension, characterized by the 
parameter R^/M^,. By fitting the limb-darkening law to a star's observed surface-brightness 
distribution, the coefficients determine the star's R-^/M^. Determining the star's i?^ using 
opt ical interferometr y and paralla x measurements gives the star's mass. The approach taken 
by iNeilson fc Lesterl (120111 . 120121 ) is quite different from the analysis done here, but the 
conclusion is the same; it is possible to determine the mass of a cool, luminous star by 
analyzing its surface brightness distribution. 

A potential interesting variation of the methods explored here is to apply them to 
radially puls ating stars, such as Cep heids. Even though Cepheids lose mass over time via a 
stellar wind (INeilson &: Lesterl |2008| ) . at a given epoch the Cepheid's mass is constant. By 
measuring the periodic variation of the star's luminosity and radius, the diagnostics explored 
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here would be able to determine the Cepheid's mass at its current evolutionary stage. 



Conclusions 



We have used the SAtlas_ODF program (ILester fc Neilsonll2008l ) to compute spherical 
model stellar atmospheres characterized by the three fundamental parameters L^, and 
i?*, which were chosen to represent typical red giants and supergiants. Because it is possible 
to observe and R^, we held the luminosity and radius constant and varied the mass. 
Searching the radiation of the models for features that varied with M^, we found that there 
is a systematic change in the surface-brightness distribution, which we characterize by the 
renormalized fractional radius, r/i^umb, at the stellar disk's half-intensity point, J(r)//(r = 
0) = 0.5. For the giant and supergiant models used here, the r/Rnrah varies by a few percent 
for the mass range we explored, which is currently below observational detection. However, 
with a larger mass variation, and a corresponding larger atmospheric extension, and with 
improving observational techniques, this signature will become viable. 

We also synthesized the flux spectra of the red giant and supergiant spherical models 
with a spectral resolving power of 10^, taking into account the spherical radiative transfer. 
The most interesting changes with mass were found to be the lines of [O I], which become 
weaker with increasing mass, and numerous lines of MgH in the X A^U band head 

around 500 nm, which strengthen with increasing mass. Both of these spectral features have 
been used previously as indicators of log^Q^g), and they can be used now to determine stellar 
mass with current spectral capabilities. 

This work has been supported by a Discovery grant from the Natural Sciences and 
Engineering Research Council of Canada to JBL, and HN acknowledges funding from the 
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explain them more clearly. 
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